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Abstract 

We address the controversial hot question concerning the validity of the loose coupling versus the lever-arm 
theories in the actomyosin dynamics by re-interpreting and extending the phenomenological washboard potential 
model proposed by some of us in a previous paper. In this new model a Brownian motion harnessing thermal 
energy is assumed to co-exist with the deterministic swing of the lever-arm, to yield an excellent fit of the set of 
data obtained by some of us on the sliding of Myosin II heads on immobilized actin filaments under various load 
conditions. Our theoretical arguments are complemented by accurate numerical simulations, and the robustness 
of the model is tested via different choices of parameters and potential profiles. 
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1 Introduction 

In 1985 an experimental set-up to measure the distance traveled by an actin filament interacting with a Myosin II 
filament 1 during a complete ATP cycle was described [1]. Under low load conditions on the actin filament, trav- 
eled distances of up to 60 nm were observed. This value was not understandable according to the most popular 
idea of the sliding mechanism, rotation or tilting of the myosin head bound on actin molecule, conceived as di- 
rectly coupled with the ATP hydrolysis (see [2]). Independent, successive experiments (see [3], [4] and references 
therein) then showed distances ranging between 4 and 10 nanometers. Such smaller range of traveled distances is 
consistent with the lever-arm theory (see [5] and [6]) that is still widely believed to account for the generation of 
the force responsible for the actin filament sliding. The level-arm theory is strictly deterministic: Each ATP cycle 
generates one single power-stroke that causes a sliding of constant length and preassigned direction of the actin 
filament. Within such framework, a tight coupling between ATP cycles and protein movements is envisaged. 

*Work performed within a joint cooperation agreement between Japan Science and Technology Corporation (JST) and Universita di Napoli 
Federico II, under partial support by MIUR (PRIN 2003) and by Campania Region. 
'Hereafter we shall use the word "myosin" as an abbreviation of "Myosin II". 
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Even though the apparent disagreement of the above mentioned difference of travelled distances could be 
understood on the base of a low duty-ratio characterizing Myosin II heads (see [8], pages 221-226), successive 
experiments by Kitamura et al. ([9]) raised anew the question of such a disagreement. Indeed, striving towards 
increasingly accurate measurements of the traveled distances during an ATP cycle, it was possible to achieve 
great improvements of the technological set-up that included design and construction of "home made" highly 
sophisticated devices. By exploiting such a technology, the experimental setup was as follows: A single myosin 
head was attached to the tip of a glass microneedle and placed near an actin filament that had been previously 
immobilized on a microscopy slide by means of optical tweezers. The deflections of the needle with respect to its 
resting position were then measured and recorded. From the obtained traces, the following three features emerged 
that are in evident disagreement with the tight coupling theory: 

1. The total traveled distance (i.e. the total displacement) of the myosin head is not constant and it can be as 
large as 30 nm. 

2. This displacement is the sum of a random number of single "steps", the amplitude of each of which equals 
the distance (~ 5.3 nm) between two successive actin monomers. During the time elapsing between two 
successive steps, the myosin head randomly jitters around an equilibrium position. 

3. Steps mainly occur in a fixed "forward" direction, although some of them occasionally take place in the 
opposite "backward" direction. Hereafter, forward steps will be taken as positive and backward steps as 
negative. The total displacement is thus the algebraic sum of the number of performed forward and backward 
steps. 

Such evidence contrasts the one-to-one relation hypothesized between ATP hydrolysis and the occurrence of the 
mechanical event consisting of the power stroke in the myosin head. Furthermore the observation of the existence 
of random elements leads one to conjecture that a significant role could be played in such context by the thermal 
agitation of the environmental molecules of the watery solution in which the involved proteins are embedded. 
This is the motivation for the assumption of the existence of a loose coupling between ATP cycle and actomyosin 
dynamics (See, for instance, [2]). 

While referring to [ 1 0] for a lucid outline of the origin of the controversy existing among the supporters of the 
loose coupling approach and the community of those faithful to the lever-arm theory, and therefore also to the tight 
coupling vision, in the present paper we shall extend the phenomenological model earlier proposed by some of us 
([1 1]) by including in it the effect of the swing of the lever-arm. Thus doing we show that the experimental data of 
myosin sliding obtained in [9] under various load conditions can be very well accounted for. 

Specifically, via the comparison of our theoretical results with those yielded by our experiments, we shall 
test the assumption that during the time interval elapsing between the attachment of myosin head to the actin 
filament and the final release of the phosphoric radical (rising phase in the sequel), the position of the myosin head 
is determined, both by the lever-arm swing and by the action of Brownian motion, that includes a macroscopic 
deterministic force responsible for a non-zero average net displacement. Such a direction-orienting force can be 
viewed as the originated by changes of chemical states that occur within the myosin head and that are fuelled by 
energy supplied by myosin head itself. Such a view is coherent with the notion of "effective driving potential" in 
the sense, for instance, of Wang and Oster [12]. As a consequence, the coupling between the ATP cycle and the 
mechanical effects should appear to be less rigid, and the actomyosin dynamics should definitely exhibit variability 
features of the kind mentioned in the above l.-j-3. items. 

The above outlined conjectured random dynamics will be described in details in Section 2 with a specific 
reference to an earlier paper [11]. Here we limit ourselves to showing that, within such a framework, it is possible 
to handle a key point of the mentioned controversy. Indeed, in [ 1 0] the essential relevance of the role played by 
the length of the myosin neck is stressed, since the lever-arm theory makes the sliding distance less than, and 
somewhat proportional to, such length. This would for instance imply that reducing the length of the myosin neck 
to a half should reduce to a half the sliding distance, which appears to be confirmed by the experiments in [ 1 3] . On 
the contrary, certain performed experiments only show slight changes of the overall displacement of the myosin 
head even after complete removal of its neck [10]. 

In order to reconcile these evident discrepancies, we shall assume that the displacement X of the myosin head 
during the ATP cycle could in general be envisaged as a linear combination as follows: 



X = rX R + dX D 



(1) 
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where Xr denotes the displacement induced by the biased thermal effects and Xu the displacement generated by 
the power stroke, and where r and d are constants, each of which hereafter will be taken as equal to or 1 . Note 
that the case r = and d = 1 yields the lever-arm theory; instead, setting r = 1 and d = depicts a purely 
random situation, in the absence of any sliding due to the power stroke. Finally, the case r = d = 1 leads to an 
integration of the two theories. With the choice r = d = 1, the controversial results related to the length of the 
neck of myosin head can be overcome by assuming that the random displacement is a few times larger than the 
deterministic one. This assumption implies that only slight changes of the total distance traveled by myosin head 
would be observable when the contribution Xjy due the power stroke is made somewhat smaller by shortening the 
length of the myosin neck. 

2 The model 

Let L denote the distance between each pair of neighboring actin monomers. As suggested in recent literature 
(see [14] and [15]) in our computations we shall take L = 5.5 nm. In addition, we shall assumed that the 
magnitude of the sliding induced by the swing of the lever-arm equals that of a step, and thus take Xjj ~ L, 

Our conjectured relation (1) with r = d = 1 is preliminarily supported by the bar chart in Fig. 4c) of [9] 
showing that at least one step is performed by each and every myosin head during the rising phase. The second 
column of Table 1 shows the heights of the columns of the mentioned bar chart, whereas the first column indicates 
net step numbers, i.e. the integral part of the ratios X/L. The quoted experiment was performed under low load 
conditions by means of microneedles having stiffness less than 0.1 pN/nm. 

Table 1 : Observed distribution of net number of steps performed by myosin heads as 
indicated in [9]. Here [X/L\ is the total step number during the entire rising phase. 



[X/L\ 


Observed frequency 




Theoretical frequency 


1 


14 





15 


2 


21 


1 


22 


3 


18 


2 


17 


4 


10 


3 


8 


5 


3 


4 


3 


6 





5 


1 


total 


66 


total 


66 



As shown by the first two columns of Table 1, one sees for instance that out of 66 observed myosin heads 
(all attached to a glass microneedle of stiffness less than 0.1 pN/nm) 14 performed a unit net steps, 21 a net 
step number equal to 2, etc., to conclude that 3 of them have performed a net step number equal to 5 implying 
total displacement of about 30 nm. Columns 3 and 4 of Table 1 show that \_Xr/L\ = [X/L\ — 1 is well fitted 
by a Poisson distribution with parameter r) = 1.5 given by the ratio of the total number of performed net steps 
(14 - + 21-1 + .. . + 3- 4 = 99) to the number of considered myosin heads (66). 

The agreement between experimentally observed frequencies and those predicted via the Poisson distribution, 
jointly with the rarity of the backward steps, leads one to conclude that the "dwell time" (namely the time interval 
elapsing between two successive steps) should be, to a good approximation, exponentially distributed. Such con- 
clusion is experimentally supported by the data summarized in Fig. 4b) of [9] leading to the estimate of about 5 ms 
for the mean dwell time. 

Our facing sequences of rare events with exponentially distributed interarrival times is strongly suggestive of a 
first-exit problem out of an interval for continuous Markov processes possessing an equilibrium point sufficiently 
far from at least one of the end points of the diffusion interval (See, for instance, [16]). 

On the ground of all foregoing considerations, with reference to the random part of the rising phase 2 we are 
led to construct a model for the actomyosin dynamics that is based on the following assumptions: 

2 By "random part of the rising phase" we denote the time interval during which the myosin head is subject exclusively to Brownian motion. 
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(i) The complex M.ADP.Pi + Energy is viewed as a point-size particle moving along an axis X on which the 
abscissa x denotes the displacement of the particle from the starting position. The positive direction of X is 
that of the forward steps of myosin head. 

(ii) The particle is embedded in a fluid. Hence, it is not only subject to a dissipative viscous force character- 
ized by a drag coefficient /?, but also to microscopic forces originating from the thermal motion of the fluid 
molecules. On account of the fluctuation-dissipation theorem, such microscopic forces can be macroscopi- 
cally described by means of a Gaussian white noise having intensity 2f3k B T, where k B is Boltzmann constant 
and T the absolute temperature. 

(iii) The global interaction of the particle with the actin filament is synthesized in a conservative unique force 
deriving from a potential U (x). The structure of the actin filament suggests that U (x) be a periodic function 
with period equal to the distance L between pairs of consecutive actin monomers: 

U{x) = U(x + rL), VreZ. (2) 

Henceforth we shall denote by La (0 < La < L) the minimum of U{x), assume U(La) = and denote by 
Uo := U (0) = U(L) the depth of the potential well, namely the maximum of U (x). 

(iv) The particle's dynamics is described by Newton's equation in which the total acting force is the sum of two 
terms: The first term is a deterministic force generated by the potential U (x) and by the viscous force, while 
the second term is the random force due to the presence of the Gaussian white noise. 

(v) Two more (constant) forces, denoted by Fi and F e , act on the particle. We assume that Fi is generated by a 
process that finds its origin in a part of the energy possessed by the particle, and take FiL <C Uq. Instead, 
F e is an external force, conceivable applied from the outside by the experimenter. 

Summing up, we are assuming that the complex M.ADP.Pi + Energy can be looked at as a Brownian particle 
subject to a tilted potential V(x): 

V(x) = U{x) - Fx (3) 

where U(x) has been indicated in (2) and F = Fi — F e . In the experimental conditions the height of the potential 
wells is Uo < 100 pN-nm, the period of the potential is L = 5.5 nm, the particle mass is to = 2.2 • 10~ 22 kg, 
the drag coefficient is j3 ~ 90 pN ns/nm and the environmental temperature is T = 293 K. Therefore, Reynolds's 
number is much less than 1 (see, also, [17]), so that the inertial term of the equation of motion can be disre- 
garded. In conclusion, the overdamped equation describing the movement of the particle is the following Langevin 
equation: 

where A(t) is a zero-mean white Gaussian noise with unit intensity, ( ' ) denotes time derivative and ( ' ) space 
derivative. 

Within this framework, this idealized Brownian particle randomly moves around an equilibrium point located at 
a minimum La of U (x). Whenever it exits the current potential well, we conventionally say that the corresponding 
myosin head has made a step, in the forward or in the backward direction according to where the exit has taken 
place. Hence, we are facing a first-exit problem of the Brownian particle from the endpoint of the current potential 
well. Taking into account the above assumptions, we conclude that the distribution of the first-exit time of the 
Brownian particle from the potential well is approximately exponential, since 

- the process described by equations (2), (3) and (4) possesses an equilibrium point at the minimum La of 

U(x); 

- the time for the particle to travel the distance La in the presence of the only force due to U (x) is (see for 
instance, [18]) r = (3L 2 A /U Q ; 

- the standard deviation of the Gaussian steady-state distribution of the process modeling the particle's motion 

is a = ^(fc B T//3) • t/2 = L A /V2^, where we have set uq = Uq/ k B T; 
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- the ratio I a '■= La/<t = y/2v>o falls well within the interval (2, 4^2) for all choices of uo € [2, 16] to which 
we shall refer in the foregoing, so that the exponential approximation for the first-exit time is valid [16]. 

In the sequel, we shall exploit the known formulas [ ] for the conditional probability p that the particle moves 
a step forward after leaving the current potential well, 

P = 1 + cxp {-FL/KT) ' (5) 
and for the mean first-exit time fx from a potential well: 



3 Determination of Fi and Uq 

Hereafter, we shall view as constants the parameters (3 and T that characterize the thermal bath and the period 
L of the potential U(x), their values being given in Section 2. Hence, the quantitative specification of the model 
described by Eqs. (2), (3) and (4) requires that numerical values be attributed to three more parameters: the depth 
Uo of the potential well, the position La of the minimum of U(x) in (0, L) and the internal force Fi. The numerical 
specification of these parameters can be performed after the function U (x) has been chosen. We shall preliminarily 
take U(x) as a symmetric (La = L/2) sawtooth potential (see Table 6). Successively, we shall test the robustness 
of our model by assuming alternative potential functions, henceforth called "potential profiles". It should be noted 
that Fi is somewhat related to the largest force that myosin is able to endogenously generate. Here, we shall not 
attempt to provide any biological justification of its genesis. We limit ourselves to pointing out that elsewhere [20], 
where similar experiments were performed on Myosin VI, it is conjectured that the week binding between actin and 
myosin is a source of distortion of the geometry of the two helixes in the actin filament. Such a distortion exposes 
the hydrophobic region of actin to myosin head, thus generating a tilt of the potential, and hence a constant force 
Fi, The existence of a tilt of the potential, conjectured in [11], has been successively supported by means of 
simulations in [21] where it is shown that in the absence of such a tilt experimental available evidence on the 
myosin motion in the presence of contrasting applied loads cannot be accounted for by any of the other models 
therein considered. Within our strictly phenomenological framework the existence of this force Fi is supported 
by Eq. (5) showing that in the absence of external applied forces (i.e. F e = 0) steps on either directions would 
be equally likely unless Fi ^ 0, in contrast with the experimental evidence on the high degree of directionality 
exhibited by motion of the myosin head. 

To proceed along the quantitative specification of our model in a way to be able to attempt the fitting of available 
experimental data, the values of Fi and of Uq must be specified. This will be done by making use of the available 
experimental data [22] shown in Table 2 and in Table 3 3 . From them it is evident that the frequency p of forward 
steps decreases as the applied load increases, while the dwell times increase with the load. 

For some fixed values of internal force Fi, use of Equation (5) has been made to calculate the theoretical 
probabilities of the particle's exit from the current potential well to the next well (i.e. the analogue of the forward 
step probabilities) as function of F e . The results are shown in Figure 1, where eight realistic values of Fi have 
been chosen in the interval [l.OO pN, 1.90 pN] . Vertical lines indicate the 3 load intervals of Table 2, whereas 
horizontal lines indicate the 3 corresponding recorded frequencies. We see that for internal forces Fi = 1.00 pN 
and Fi = 1.90 pN the plotted curves do not meet the requirement of leading to the experimentally recorded 
frequency p = 0.8571, whatever value F e is chosen within interval [0 pN, 0.50 pN] . Similarly, for Fi = 1.55 pN 
no value of the computed probability equals the frequency p = 0.8163 for F e ranging in [0.50 pN, 1.00 pN]. 
Instead, all remaining 5 curves referring to values of Fi ranging from 1.60 pN to 1.80 pN in steps of magnitude 
0.05 pN, from below upward, are in agreement with the experimental values of Table 2. Hence, the interval of 
values for Fj to be selected in order to secure the fitting of the experimental data is [l.60 pN, 1.80 pN] . 

We now come to the estimation of the last parameter, Uq, i.e. of the depth of the potential well. This is done 
by exploiting Eq. 6 in which the left-hand side is viewed as the function \x = /u(Fj, F e , Uq; L, (3, T) and fixed to 

3 Note that loads and dwell times in Table 2 and 3 must be viewed as averages to which confidence intervals are associated. For instance, 
the load 0.046 pN is the result of all measurements for which the product of the stiffness of the glass microneedle times the distance traveled 
by the myosin head falls around 0.046 pN. The corresponding dwell time 5.3 ms is to be viewed as the arithmetic average of the dwell times 
recorded during these measurements. 



4 NET STEP NUMBER DISTRIB UTION 



6 



Table 2: For three conditions of the applied load Cmeasured at 
the end of the rising phases, the table lists the recorded numbers 
fif of forward steps, fib of backward steps, the total number of 
steps and the percentage p of forward steps. 



Table 3: Recorded dwell times p. for 
different values of the applied load 
C. 



G (pN) 


% 


hb 


fif + m 


h f 

P= - . . 
Uf + n b 


] 0.0, 0.5] 


54 


9 


63 


0.8571 


] 0.5, 1.0] 


40 


9 


49 


0.8163 


] 1.0, 2.0] 


29 


19 


48 


0.6042 



C (pN) 


u (ms) 


0.046 


5.3 


0.190 


5.7 


0.300 


6.0 


0.470 


7.1 


0.690 


8.9 


0.830 


6.2 


1.240 


11.1 


1.890 


11.0 



5.3 ms, that (see Table 3) corresponds to the smallest applied load 0.046 pN. Since temperature T, period L of 
the potential V(x) = U (x) — (Fi — F e )x, and drag coefficient (3 are specified, if we take F e = 0.046 pN Eq. (6) 
makes Uq an implicit function of Fi. 

Note that as Fi increases (i.e. as the potential tilt increases) the mean first-exit time fj, decreases. Hence, in 
order to keep fi constantly equal to 5.3 ns, while Fi ranges in the interval [1.60 pN, 1.80 pN] , depth Uq must be 
taken as a monotonically increasing function of Fi, From /i(1.60, 0.046, Uq) — 5.3 and ^(1.80, 0.046, Uq) = 5.3, 
Eq. (6) yields Uq w 15.632 k B T and U w 15.755 fc B T, respectively. 

In order to test the agreement of our model with the experimental values of mean dwell times for the various 
loads (see Table 3), we make use of Eq. (6) to determine fi as a function of F e for the pairs (1.60 pN,15.632 k B T) 
and (1.80 pN,15.755 fc B T), involving the extrema of the determined values for Fi and Uq. The obtained values 
are showing Figure 2, where the corresponding experimental mean dwell times are also indicated. Obviously, all 
other pairs of admissible values of Fi and Uq lead to curves lying inside of the above two pairs. Inspection of 
Figure 2, jointly with the magnitudes of the confidence intervals associated to each experimental value [15], shows 
that for small values of F e the agreement between experimental data and theoretical predictions is good. Indeed, 
up to the first 3 values of Table 2 the agreement is excellent, to become more than acceptable for larger values 
of F e up to 0.83 pN. The discrepancies shown by the remaining two experimental data are a consequence of the 
crude assumption (vi) above by which the applied load is assumed to be strictly parallel to the direction of motion. 
This is acceptable for reasonably small loads, but certainly unrealistic for large loads: In this case a significant 
orthogonal component should be expected, whose effect is alike to an increase of the depth of the potential well 
Uq. Nevertheless, in the interval between the last two load values of Table 2, a qualitative similar behavior of 
theoretical curves and experimental values is present. 



4 Net step number distribution 

In the present Section we shall implement the model described by Eqs.(2)-(4) in order to obtain, via a suitable 
simulation procedure, the distribution of the net number of steps performed by the myosin head during the time 
interval elapsing between the instant when the ATP molecule is hydrolyzed and the final release of the phosphoric 
radical, i.e. during the random part of the rising phase. The results of our simulations will then be discussed with 
reference to the experimental data of Table 1 . We recall that such data refer to 66 myosin heads that have globally 
performed 99 net steps. 

Our simulation procedure is based on a discretized version of Eq. (4) and on a routine for generating Gaussian 
pseudorandom numbers in a way to determine step by step the positions achieved by 66 Brownian particles all 
originating at xq at time 0. The traveled distances at the end of the individual random rising phases are recorded, 
which finally leads one to the determination of the 66 net step numbers. Such procedure has been implemented on 
an IBM SP4 parallel supercomputer and repeated 1600 times to obtain a reliable statistics. 

In order to solve Eq. (4) numerically, one has to specify preliminarily initial position xq of the Brownian 
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0.0 0.5 1.0 1.5 2.0 

External applied force (pN) 

Figure 1 : The conditional probability p that the particle moves a step forward after leaving 
the current potential well is plotted as function of the external applied force F e for various 
values of internal force F, (pN). Temperature and potential period are taken as T = 293 K 
and L = 5.5 nm, respectively. The three horizontal lines whose ordinates are 0.8571, 0.8163 
and 0.6042, respectively, denote the sample percentages p of forward steps under the three 
load conditions indicated in Table 2. 

particle and duration of each sample path. We shall safely take xq = L/2 since, whatever the actual initial 
position of the particle inside of the potential well, the its relaxation time is much smaller than the mean first-exit 
time. Incidentally, we note that such a choice is consistent with the procedure adopted in Kitamura et al. [9] (see 
legend to Figure 2d), where to determine the average recorded trajectory, the starting positions of the rising phases 
have been synchronized. The specification of sample path duration is somewhat more involved because is a 
random variable of unknown distribution. Nevertheless, we can appeal to the approximate exponential distribution 
of the first-exit time, motivated by the depth of the potential well, to assume that is approximately gamma- 
distributed, with probability density T(n, v) where [i is the mean first-exit time from a potential well and v is the 
mean total exits of the Brownian particle. While /i is obtained via Eq. (6), our estimate v of v is obtained by using 
the data of Table 1 and Eq. (5): 

99 

V ~ 66(2p- 1)' 

Finally, for each specified potential well Uq, the size At of the time parsing in Eq. (6) is determined by pro- 
gressively reducing it until the obtained distribution becomes appreciably invariant. For an immediate comparison, 
Table 4 shows the distribution obtained via simulation for a potential Uq = 5 k B T and a parsing time At = 0.25 ns 
jointly with the experimental distribution. The agreement between observed and numerical frequencies is more 
than satisfactory. The unique somewhat apparently appreciable discrepancy is that 2 (out of 66) simulated Brown- 
ian particles are seen to perform 1 negative net step number, implying a zero total displacement due to the super- 
position of the effect of the final power stroke. However, at the present stage, such a discrepancy should not be 
taken as particularly significant. Indeed, the presently available experimental setup does not allow one to monitor 
simultaneously the trajectory of the myosin head and the occurrence of hydrolysis of an ATP molecule. Therefore, 
observing at least one positive net step number of myosin head indicates that one ATP hydrolysis has occurred, 
and consequently such a case has been included among the recorded data. Instead, in the case when a zero net step 
number is observed, it is impossible to claim that an ATP molecule has been hydrolyzed, and hence such cases 
have been disregarded. Future endeavours could aim at the design of an experimental setup able to prove that a 
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Figure 2: Mean first-exit time fj, is plotted as a function of the external applied force F e . 
The chosen values of Fi (pN) and Uq (k B T) are indicated for each curve. Dots represent the 
experimental dwell times of Table 2. Drag coefficient j3 has been taken as 90 pNns/nm, 
whereas temperature and period of the potential are the same as in Figure 1 . 

zero net step number can actually occur as a consequence of the hydrolysis of an ATP molecule. 

The chosen value of Uq is motivated by a twofold consideration. First of all, it is large enough to imply that 
the first-exit time distribution is approximately exponential. This is indeed supported by the numerically evaluated 
coefficient of variation that has been found to be 0.96, close to the value of 1.0 for the exponential distribution. 
Second, and most relevant consideration, is that it is reasonable to conceive that, under the assumed overdamped 
regime, the net step number is insensitive to the depth of the potential well. Indeed, the forward exit probability p 
given by (5), under the fixed environmental temperature, only depends on the energy FL, namely on the difference 
of potential V(x) over one period, thus being independent of the depth Uq. Table 5 evidently supports such 
conclusion. Indeed, it indicates that, for instance, by doubling the depth of the potential well Uq, the net step 
frequencies are not affected significantly, even for different choices of internal forces 4 . 

Note that the parsing parameter At must be determined with specific reference to the magnitude of Uq. Indeed, 
At must be much smaller than the relaxation time of the force —U'(x) oc (3L 2 /Uq. Furthermore, it must be such 
as to cope with the random forces due to the thermal bath. Therefore, the mean square displacement per unit time 
should remain constant as the depth Uq of U(x) is made to change, which implies an inverse dependence of the 
magnitude of At on the square of magnitude of the potential well. This is the motivation for the choices of At 
indicated in Table 5. 

5 The role of potential forms and asymmetries 

In [1 1] the idea and role of a washboard-potential played in myosin dynamics was introduced and exploited with 
the only reference to a parabolic potential (parabolic "profile", in the currently adopted terminology), though 
without any consideration on model robustness. This task is now accomplished in the more general framework of 
the present model, by considering not only parabolic, but also other types of profiles: sawtooth, cosine-like and 

4 The presence of non-integer numbers in Table 4, 5 and 8 is a consequence of our adopted estimation procedure. The half-width of the 
related 95%-confidence intervals has been seen never to exceed 0.2. For comparison reasons we have not rounded out the raw numbers to the 
nearest integers. 
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Table 4: Net step frequency distributions for Uq — 5 k B T and for three values of Fi chosen in the 
admissible interval. Here F e — 0, xq — L/2 and At = 0.05 ns. Other parameters have been chosen 
as in Figure 2. The last 5 rows show the coefficient of variation, probability of a step forward and 
mean exit time, all obtained via simulations. The theoretical values of p and fj,, obtained via Eqs. (5) 
and (6), are also indicated. 



[Xr/l\ 


Observed frequency 


Numerical frequency 






ti = 1.70 pN 


Fi = 1.75 pN 


Ft = 1.80 pN 


-1 




2.05 


1.90 


1.76 





14 


13.79 


13.63 


13.88 


1 


21 


20.72 


20.83 


20.93 


2 


18 


16.41 


16.43 


16.47 


3 


10 


8.53 


8.59 


8.54 


4 


3 


3.14 


3.26 


3.16 


5 





0.93 


0.93 


0.88 


CV 


0.966 


0.964 


0.963 


P 


0.909 


0.914 


0.920 


P 


0.910 


0.915 


0.920 


~p (ns) 


1235.057 


1207.288 


1180.852 


//(ns) 


1233.784 


1206.048 


1178.695 



Lindner-type (see Table 6). The effects of the natural asymmetry present in the acto-myosin system is investigated 
within our model with asymmetric sawtooth potential. A plot of the considered potentials over one period are 
shown in Figure 3, whereas Figure 4 refers to the corresponding generated forces. Lindner-type potential has been 
indicated for two values of parameter S. It is not difficult to see that S — > yields the cosine profile, whereas the 
potential flattens down in the middle as 6 increases. 

The independence of the exit probability p of the potential's profile implies that the net step distribution is 
profile-independent as well. This is also evident from Table 8, where the net step distributions are reported for 
each of the above-considered four potential profiles. 

Next task is to pinpoint the effects of the potential profiles on the mean first-exit time. To this purpose, we refer 
to Table 3 showing that the mean dwell time in lowest load condition is fj, = 5.3 ms. After choosing Fi = 1.75 pN, 
we then make use of Eq.(6) for each and every one of the four considered potentials imposing that the left hand 
side equals such value of yu. By iterated numerical integrations, for each potential profile the corresponding value 
of Uq is finally determined. The result are listed in Table 7, where the effect of the parameter 6 in Lindner-type 
profile has been detailed. 



Table 5: Net step frequency distribution numerically obtained. Here F e = and xq = L/2. Other parameters are 
chosen as in Figure 2. Parsing steps are chosen as follows: At = 0.1 ns for Uq = 4 k B T, and At = 0.025 ns for 
Uq = 8 k B T. The indicated values of Fi belong to the admissible interval. 





Fi = 1.70 pN 


Fi = 1.75 pN 


Fi = 1.80 pN 


[X R /L\ 


U = 4 k B T 


U = 8 k B T 


U = 4 k B T 


U = 8 k B T 


U = 4 k B T 


U = 8 k B T 


-1 


2.11 


2.05 


1.94 


1.90 


1.85 


1.81 





13.56 


13.96 


13.54 


13.76 


13.48 


14.00 


1 


21.13 


20.26 


20.92 


20.26 


21.30 


20.46 


2 


16.70 


16.11 


16.93 


16.29 


16.98 


16.09 


3 


8.43 


8.58 


8.53 


8.61 


8.40 


8.56 


4 


2.95 


3.40 


3.01 


3.53 


2.91 


3.49 


5 


0.75 


1.12 


0.78 


1.15 


0.76 


1.11 
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Table 6: Potentials' profiles. 



Sawtooth 



U s (x) 



-Uo 

L A 
U 



L-L 



(x-L A ), 0<x<L A 
(x — La), La < x < L 



Parabolic 



Up(x) = 



L 2 /4 



L 

x 

2 



Cosine 



U c (x) 



Uo 
2 



2tt 
[ — x 



1 



Lindner 



U L (x) = 



Uo 

>28 _ 



, 2tt 
cos I —X ) + 1 



- 1 



Table 7: For each potential profile 
the depth Uo of the potential well 
is indicated. Here Fj = 1.75 pN, 
F e = 0.046 pN and the other pa- 
rameters are the same as in Fig- 
ure 2. 



Type 


Up (k B T) 


Sawtooth 


15.723 


Parabolic 


15.043 


Cosine 


13.944 


Lindner (S = 0.1) 


13.942 


Lindner (6 = 0.5) 


13.918 


Lindner (S = 1) 


13.851 


Lindner (S = 2) 


13.697 


Lindner (S = 3) 


13.611 


Lindner {5 = 4) 


13.591 


Lindner (5 = 5) 


13.604 


Lindner {5 = 10) 


13.749 



The behavior of mean first-exit time fx as a function of the magnitude of the external force is shown in Figure 5. 
All considered cases of Table 7 lead to graphs falling in the region bounded by the lowest and the highest curves. 
We point out that changing the potential profiles never yields mean first-exit time changes exceeding 10%. Hence, 
we are led to conclude that the depth Uo of the potential well can be tuned to acceptable biological values by a 
suitable selection of the potential profile, without affecting appreciably the value of the mean first-exit time. For 
instance, switching from sawtooth to Lindner-type potential with 6 = 2, lowers Uo by more than 2 k B T. (See 
Table 7). It is thus conceivable that a variety of potential profiles exists such that the height of the potential well 
can be further lowered, in a way to switch from about 15 k B T as indicated in [21] to about 5 k B T as suggested 
in [ 1 8] . The quantitative analysis and comparison with available data has been performed under the assumption of 
rigorous symmetries exhibited by the profiles of the potentials generating the periodic force acting on the Brownian 
particles. It is, however, presumable that the complex biological reality underlying the observed motion of the 
myosin heads may require to relax the rigorous symmetry assumption. To test the effect of symmetry breaking, 
we take into consideration the sawtooth potential Us(x) of Table 6 and make it asymmetric by taking La ^ L/2, 
namely by shifting in either direction the point of minimum of the potential. Thus doing, the introduced asymmetry 
should affect the motion of the Brownian particles. While the probabilities of exit from the current potential well 
are insensitive to the potential's profile, and hence also to its asymmetry, the mean first-exit time fi is clearly 
affected by it, as shown by Eq. (6). The quantitative dependence of /j, on the potential's asymmetry is indicated in 
Figure 6, where on the abscissa the external force F e is indicated. Each curve is characterized by two parameters: 
the point La of minimum and the depth Uq of the potential. From top to bottom, the first curve is characterized 
by asymmetry La = 2.0 nm, implying a somewhat steeper rise of the potential in the backward direction. Such 
asymmetry is mitigated in the next curve {La = 2.5 nm). The third curve, indicated in bold, is used as a comparison 
tool, as it refer to symmetric profile {La = L/2). The remaining two curves are labeled by La = 3.0 nm and 
La = 3.5 nm, thus representing mirror-image situation with respect to the zero-asymmetry case. Now the forward 
edges of the potentials are steeper. For each curve, the indicated value of Uo has been determined by imposing that 
/i = 5.3 ms when F e = 0.046 pN (lowest load), so that all curves originate at the same point. As Figure 6 shows, 
changes of the potential's asymmetry of 30% in either direction do not affect greatly the magnitude of the mean 
first-exit time. Finally, we remark that suitable choices of backward asymmetry (such La = 2.0 nm in Figure 6) 
improve the fitting of the experimental data even for larger load values. 
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Figure 3: Plots of the potentials as function of the Figure 4: Conservative forces originated by the po- 

position within each pair of consecutive monomers. tentials of Figure 3. 

6 Concluding remarks 

The model of actomyosin dynamics discussed in the foregoing rests on the assumption that the total energy made 
available to the myosin head by the ATP molecule hydrolysis and by the thermal bath has a two-fold overall 
role: To produce the power stroke predicted by the lever-arm model and also to generate the kind of sliding of 
myosin head on the actin filament in the "loose coupling" mechanism originally hypothesized in [23] and then 
experimentally demonstrated in [9]. This is expressed mathematically via the representation of the displacement 
of a particle consisting of a combination of a deterministic part and of a random component. The latter is generated 
by the simultaneous presence of a washboard-type potential and a random force arising from thermal fluctuations. 
Our model has then been tested by making use of a set of data on the dwell times and step frequencies of myosin 
heads under various load conditions. We have shown that by a suitable tuning of the internal force and depth of the 
potential well, the theoretically calculated probability p and mean first-exit time /i of the representative Brownian 
particle are in good agreement with their biological counterpart. A second set of experimental data concerning 
the net step number distribution of myosin heads under low load conditions has then been exploited to show that 
the washboard potential used by us is able to reproduce such distribution within the mathematical analogy. To the 
rather small size of the experimental sample should be ascribed the 3% discrepancy represented by the two net 
backward displacements predicted by our model. 

Next, the robustness of our model has been tested by inserting 4 different potential profiles in the Langevin 
equation of motion. The consequently performed calculations have shown that the mean exit-time depends on the 
external force in a way that is essentially insensitive to the chosen potential's profile. The chosen profile, instead, 
has been seen to play an essential role in that it critically relates the depth of the potential well to magnitude of the 
mean exit-time. A finer tuning of the mean exit-time is finally achieved by regulating the level of asymmetry of 
the potential's profile. 
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Figure 5: For each potential profile the mean first-exit time /i is plotted as a function of the 
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We have taken Fi = 1.75 pN, while the other parameters have been chosen as in Figure 2. 
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Figure 6: For the sawtooth potential profile with the indicated value of asymmetry La, the 
mean first-exit time fx is plotted as a function of the external force F e . The corresponding 
values of Uq such that one obtains \i = 5.3 ms for F e = 0.046 pN, are also indicated on 
each curve. We have taken Fi = 1.75 pN, while the other parameters have been chosen as in 
Figure 2. 



